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An analysis of voids using cosmological A^-body simulations of cold dark matter 
(CDM) models is presented. It employs a robust statistics of voids, that was recently 

■ applied to discriminate between data from the Las Campanas Redshift Survey (LCRS) 

and different cosmological models. Here we extend the analysis to 3D and show that 
typical void sizes D in the simulated galaxy samples obey a linear scaling relation with 
the mean galaxy separation \: D = Dq + x X. It has the same slope i' as in 2D, but 
with lower absolute void sizes. The scaling relation is able to discriminate between 

OO . different cosmologies. For the best standard ACDM model, the slope of the scaling 

•/^ ' relation for voids in the dark matter halos is too steep as compared to the LCRS, with 

too small void sizes for well sampled data sets. 

By considering a range of CDM simulations we further investigate the scaling 
relation for voids within the distribution of dark matter halos and other properties of 
underdense regions. The scaling relation of voids for dark matter halos with increasing 

■ mass thresholds is even steeper than that for samples of galaxy-mass halos where we 
Oh' sparse sample the data. This shows the stronger clustering of more massive halos. 

Further, we find a correlation of the void size to its central and environmental average 
^ • density. We study the evolution of the void size distribution of dark matter halos at 

^ I redshifts up to z = 3 measuring the sizes in comoving coordinates. While there is 

Ci ■ little sign of an evolution in samples of small DM halos with Vdrc ~ 90 km/s, voids 

in halos with circular velocity over Vdrc = 200 km/s are larger at redshift z = 3 due 
to the smaller halo number density. The flow of dark matter from the underdense 
to overdense regions in an early established network of large scale structure is also 
^ ' imprinted in the evolution of the density profiles with a relative density decrease in 

void centers by A{p/p) ~ 0.18 per redshift unit between z = 3 and z = 0. 

Key words: cosmology: dark matter - galaxies: formation - large scale structure of 
the universe. 



1 INTRODUCTION 

Voids in the large scale structure of the universe are the fo- 
cus of the present paper. The building blocks of the large 
scale structure in the universe are superclusters, filaments 
of galaxies, and great walls. These structure elements are 
separated by regions where the density of visible matter is 
much less than the average density. Although several in- 
vestigations of voids have been carried out (cp. Miiller et 
al. 2000 and references therein), some key questions remain 
open. Voids were found in different sizes and properties, de- 
pending on the objects observed. The lower space density of 
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clusters and groups of galaxies leads to larger voids as com- 
pared to galajcy samples. Of particular interest is the ques- 
tion of the fraction of the total volume contained in voids. 
Obviously the answer depends not only on the completeness 
and quality of the data, but also on the definition used for 
voids. Several authors differ in the way that voids are identi- 
fied. While in early papers (e.g. Einasto et al 1983 and Oort 
1983) a qualitative description was used, later authors used 
spherical or spheroidal shaped regions which were fitted into 
the low density regions, e.g. Einasto et al. (1989), El- Ad & 
Piran (2000), Plionis & Basilakos (2001) and most recently 
Hoyle & Vogeley (2001a,b). Aikio and Mahonen (1998) pro- 
posed a different approach using points of maximal distance 
to their nearest neighboring galaxy. Voids detected with this 
algorithm are arbitrary shaped, nearly convex regions con- 
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taining no galaxies as they are formed from a collection of 
subvoids. This method yields voids which are very similar to 
the results of our algorithm described below. These methods 
have the advantage of covering almost the total volume of 
the sample with voids, so one can determine robustly the 
volume fraction that is contained in voids of different void 
sizes. Consequently the volume fraction of voids is a func- 
tion of the void size rather than a pure number as discussed 
below. 

Using cosmological A'^-body simulations we can explore 
some of the selection effects arising while investigating data 
with complicated survey geometries and object selections. 
In a simulation the structure that forms and evolves due to 
gravitational clustering is fully known and accessible. There- 
fore it is possible to study properties of voids such as size, 
abundance and density of dark matter and small galaxy ha- 
los in the voids and in the environment of the voids. A draw- 
back of this approach is that the relation of the dark matter 
distribution to visible matter such as galaxies and clusters 
of galaxies is still under investigation. A careful study of 
the underdense regions in dark matter density fields and a 
comparison with observed voids will allow us to gain further 
insight. One challenging idea that can be realised in the sim- 
ulations is to run the clock backwards and make predictions 
concerning the void properties in the past at higher red- 
shifts. The effect of a possible evolution is one of our main 
areas of interest. 

Most of our analysis is concentrated on the presently 
most popular set of the cosmological parameters, the ACDM 
model. This model includes a cosmological constant Q.a = 
0.7 with a matter content Q.m =0.3. The resulting universe 
is spatially flat, a preference of inflation theory. There are 
presently several observations which cover different aspects 
of modern cosmology and independently favor this model 
among other cold dark matter scenarios. In particular we 
cite the rece nt measurements of the light curves of distant 
supernovae ( Goldhaber et al. 2001 ), measurement of the 
acoustic peaks in the spectrum of the anisotrop ics of the 
cosmic microwave background (Jaffe et al. 2001), and the 
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We propose to undertake an additional test of the pop- 
ular ACDM model, by comparing abundance and properties 
of the voids that occur in large A''-body simulations of this 
model against existing observations. Recently in an analy- 
sis of the Las Campanas Redshift Survey (LCRS) a reliable 
method of determining the voids, together with a robust 
statistics in terms of median and quartile values of the size 
distribution, has been developed and ap plied to discriminat e 
between different cosmological models ( Miiller et al. 2000 ). 
There we studied the area fraction of voids in an unambigu- 
ous way. Voids in the previous literature are often defined as 
empty or nearly empty regions exceeding a certain arbitrary 
chosen minimum size. To be independent of such a mini- 
mum size, we identify voids of all sizes and let them occupy 
the total volume of the survey or the sample under study. A 
cumulative size distribution shows which fraction of space 
is covered by voids up to a certain size. This cumulative 
size distribution has a very similar shape in different parts 



of the LCRS despite differences in sampling and possibly 
cosmic variance. It represents the self similarity of the void 
distribution and allows the definition of median and quartile 
values of void sizes. 

Previously (Miiller et al. 2000) we showed that the me- 
dian and quartile values of the cumulative void size distri- 
bution in well defined samples obey a linear scaling relation 
with the mean separation of galaxies in the studied sample. 
The slope of this relation is an important characteristic for 
the clustering properties of the observed large scale struc- 
ture, and it proves to be able of discrimination between cos- 
mological models. The analysis of the voids in the volume- 
limited subsamples of the LCRS and the comparison to the 
mock samples constructed from different CDM models in- 
dependently favors a model with positive vacuum energy 
ACDM (and/or models with a break in the primordial power 
spectrum of density fluctuations, Kates et al. 1995). 

Here we want to extend this analysis to 3 dimensions 
and make use of the full information contained in large A^- 
body simulations. In prospect of the upcoming surveys of the 
large scale structure, testing the expectations of the models 
will be possible in the near future. One of our aims is to 
trace the evolution of underdense regions in the universe, 
and closely related, the evolution of the void size distribution 
using simulated galaxy halos. These were obtained in various 
iV-body simulations of the ACDM and the OCDM model 
with various box sizes and resolutions. Therefore we are able 
to explore the scope between large simulation boxes and high 
spatial resolution. We study the connection of the void size 
distribution to the network of high and low density regions 
of the large scale matter distribution. Thereby we extend the 
previous analytic and numeric approach of van de Weygaert 
& van Kampen (1993). Also we investigate the growth of 
the size of the underdense regions as well as the evolution 
of their density proflles. 

The present paper is organized as follows. In the next 
Section we describe the method we adopted, and in Section 
^ wepresent further analysis of voids in the LCRS. Sec- 
tion y contains a description of the CDM models and the 
simulations we use for our analysis. Section |^ is devoted to 
the study of voids and underdensity regions at present time 
2 = according to A''-body simulations. The evolution of 
voids is finally examined in Section ^ by investigating the 
cumulative void size distribution at different epochs and the 
change of the corresponding averaged density profiles of the 
voids. We close with a discussion and give some prospects 
for the future. 



2 SEARCHING AND QUANTIFYING VOIDS 

The algorithm used for the identification of voids i s anal- 
ogous to t he one employed in our previous paper ( Miiller 
et al. 200C ) . It is an adapted version of the void search algo- 
rithm originally proposed by Kauffmann and Fairall (1991). 
The principle of the algorithm can be sketched as follows. A 
density field is realised on a high resolution grid where each 
galaxy occupies one grid cell. First the largest void is found 
and the volume it occupies is assigned to it, so that no other 
void can enter this region. It proceeds to the next largest 
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void and continues until there remains no more empty space 
for voids. Since the number of occupied cells usually rep- 
resents a negligible part of the total number of cells, the 
total volume in voids occupies almost the total volume of 
the survey. 

Each time the void finder identifies a void, it starts with 
placing the largest possible base void which is a cube of 
empty cells with maximal length to the grid of the den- 
sity field. The next step of the search goes along the six 2- 
dimensional faces of the base voids. Accumulations of empty 
cells are attached to the faces of the cubic shaped base void. 
The condition for this is that the area covered by empty 
cells is larger than two third of the face area. Tests show 
that there is no significant dependence on this condition. 
This step is repeated taking each added extension as the 
parent face for the next extension as long as the area condi- 
tion holds. The procedure of adding extensions starts with 
placing the largest square of empty cells on this area and 
attach rows of empty cells to each side of the square. Also 
the length of the row has to exceed two third of the length 
of the base square. In the next step again the attached row 
is taken as the parent for the next extension along its own 
edge. After the end of this procedure, the empty cells of the 
identified void are marked, so that subsequent voids will not 
enter or connect to the region assigned to this void. 

In case of the two dimensional surveys such as the LCRS 
subsamples, the basis of the void search is a square laying in 
the plane of the data. Therefore the algorithm starts with 
the second step, i.e. only the part with the squares and the 
rows is applied here. 

Physical motivation of this procedure comes from the 
fact that underdense regions in the universe tend to evolve 
into a nearly spherical shape with time (cp. e.g. van de 
Weygaert & van Kampen 1993). While overdensity regions 
decouple from the cosmic expansion and collapse to pan- 
cakes, filaments and clusters, the underdense regions grow 
faster than the background Hubble expansion. This has the 
effect that the initial geometry, which might be arbitrary, 
evolves to a more convex shape during later epochs. The 
above scheme takes this into account in two steps. At first, 
the void search starts from cubic base voids. Then it adds 
extensions in a way that approximately conserves convex- 
ity. The algorithm identifies neighboring voids separately if 
a galaxy between them stops the extensions. The voids are 
not connected by tunnel or finger like structures. A tunnel 
between two compact voids is added to the larger void. It 
may enter the smaller void region and subdivide it, but this 
extension will not enclose the smaller void. The algorithm 
will later find the empty parts of the second smaller void 
and identify them as separate void regions. Therefore the 
present approach recognizes large voids more robustly. 

In density fields constructed from redshift surveys or 
A''-body simulations small voids are much more abundant. 
Because data quality and sample completeness can affect 
each individual void, in particular its size, shape or even 
existence, we need to define statistical values, representing 
large, medium sized and small voids. In order to study the 
void size distribution with samples of different density, we 
often dilute the simulated samples, and the stochastic di- 



lution mimics the effect of sparse sampling in the observa- 
tions and affects the occurrence of voids. As it was shown 
by Sheth (1996) the process of random dilution is linked to 
higher order correlations of the large scale galaxy distribu- 
tion, in particular in the application to counts in cells statis- 
tics. Large voids are less probable to be produced purely by 
random dilution or sparse sampling. We regard them as more 
prominent for the large scale structure geometry. 



3 2D VOID ANALYSIS OF THE LCRS 

Different from our recent approach in the analysis of the 
LCRS, we do not use the length of the base void as the 
characteristic size of the voids. Often the volume contained 
in the extensions is comparable to the base void. Therefore 
we take this into consideration and define the diameter D of 
a sphere as an effective void size, that has a volume equal to 
the total volume of the void, base volume plus extensions. 
In case of the 2D LCRS data the diameter of a circle is 
used. Further we calculate a cumulative void size distribu- 
tion, which shows how much of the total volume of the sam- 
ple is contained in voids up to a certain size. This cumulative 
distribution is a statistically robust description of the voids 
as was shown in Miiller et al. (2000). As representatives for 
typical void sizes in the galaxy distribution we compute the 
median and quartile values of the area fraction of voids for 
all samples. In Miiller et al. (2000) it was shown that these 
values obey a linear scaling relation D = Do + v x X. The 
mean galaxy separation A varies due to different sampling 
fractions in the analysed data and due to different galaxy 
densities for changing magnitude ranges. We always select 
volume limited data sets by imposing absolute magnitude 
limits. The residuum size Do - the offset of the void size at 
formally zero galaxy separation A = - describes the devi- 
ation of the scaling relation from pure proportionality. The 
galaxy separation varies due to different galaxy densities for 
the absolute magnitude cuts in selecting volume limited data 
sets, and in different sampling fractions in the LCRS slices. 

To test the robustness of the scaling relation in the 
LCRS found in Miiller et al. (2000), here we add 3 more 
analyses to the original one; 

(1) The effects of varying sampling fraction in the data. 
Each of the 6 LCRS slices consists of more than 50 observa- 
tional fields with varying fractions of measured galaxy red- 
shifts. To obtain homogeneous sampling fractions, in our 
previous paper we randomly diluted all fields to the sam- 
pling fraction of the field with the minimum sampling. Here 
we repeat this approach 10 times to evaluate the variance of 
the void size distribution. The result and the fitted scaling 
relations are shown in Fig. |l| and in Table |l|. The accuracy of 
the fits and the narrowness of the error ranges are remark- 
able. The different LCRS slices sample the voids in different 
parts of the universe homogeneously and lead to stable re- 
sults. This is the basis for our cosmological conclusions. 

(2) The effects of the grid orientation. 

In addition to the original paper we treated systematic er- 
rors introduced by the grid orientation in the determination 
of the median and quartile sizes D. These were estimated 
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Figure 1. Median (filled circles), upper and lower quartile values 
(open circles) of void sizes in the LCRS volume limited samples. 
The error bars show the variation of the results obtained by re- 
peating the selection of the subsamples. 



by repeating the void search algorithm while the volume- 
limited LCRS samples were rotated with respect to the grid 
in steps of 10°. In Fig. ^ we show the scatter in median and 
quartile values due to the grid orientation by the vertical 
error bars on the diamonds and squares of the 14 volume 
limited data sets of the LCRS. Note that these data lie on a 
narrow range around the fits to the scaling relations shown 
by solid lines. 

(3) The effect of further random dilution in homogeneous 
data sets. 

We also studied the scaling relation for subsamples that stem 
from well sampled homogeneous datasets with mean galaxy 
separation A < 12/!,~^Mpc, that were further randomly di- 
luted. Fig. 1^ shows that the median and quartile void sizes of 
these data obey the same scaling relation (dashed lines) as 
in the originally studied data sets. The scatter is somewhat 
higher than the original data, which is due to the strong 
dilution of the galaxy numbers in the restricted volumes, es- 
pecially concerning the median and the upper quartile. We 
take this as an upper bound on the variance of the void size 
distribution. These data underline that the scaling relation 
results from different galaxy densities in different samples 
and not from different galaxy properties as magnitude lim- 
its or galaxy types. The last point was discussed in detail in 
Miiller et al. (2000). 

The parameters of the scaling relations for the effec- 
tive void sizes D are given in Table |l] first for the LCRS 
data using the different random selections described in point 
(1) above to estimate the error bars. The next row in the 
Table contains the parameters of the scaling relation using 
different grid orientations. The third row gives the results 
of a linear fit to the data including the diluted subsamples 
shown in Fig. The three parameter sets agree remarkably 
well within the errors. However the scatter among the last 
sets of data is much larger and a test of the fits implies 



Figure 2. Median (large squares), upper and lower quartile val- 
ues (diamonds) of void sizes in the LCRS volume limited samples 
with rotated grid of the density field. The scaling relation for 
these data sets are shown as solid lines. The quartile values for 
diluted subsamples of rich data sets (A < 12/i~^Mpc) are given as 
open and filled circles with the scaling relations as dashed lines. 



that only the errors contained in the first two lines can be 
regarded as reliable. For a comparison with our earlier pa- 
per, we also give the parameters for the scaling relation of 
the base void lengths. The absolute values are smaller be- 
cause the base void lengths are smaller than the effective 
void sizes by about \/2. The new analysis better takes into 
account the space covered by the extensions to the base void, 
as described above in Section |^ 



4 CDM SIMULATIONS 

To evaluate the void size distribution we performed a set 
of numerical simulations of realistic CDM co smologies. As 
already shown previously (Muller et al. 2000), the void size 
distribution and the scaling relation depend on the cosmo- 
logical parameters, and in particular on the power of the 
primordial fluctuation spectrum on large scales as described 
by the shape parameter F = Mlm, where h is the dimen- 
sionless Hubble parameter h = _ffo/100 km/s/Mpc and Q.m 
the dimensionless matter density. Results that best repro- 
duce the observed void size distribution were obtained for a 
standard ACDM model and for a model with a steplike pri- 
mordial power spectrum, i.e. with a similar spectral shape 
as the widely discussed rCDM model. But even there it re- 
mained a significant discrepancy in the abundance of large 
voids and in the scaling relation of the void size distribution 
with the mean galaxy separation that becomes significant for 
the well sampled LCRS data sets. There we could not answer 
the question, whether a scale dependent bias can influence 
the void statistics, and in particular, what is the effect of 
the suppression of the galaxy formation probability in un- 
derdense regions. Such a suppression is som etimes denoted 
as large-scale bias (Doroshkevich et al. 199!:). Furthermore, 
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Table 1. Void size scaling relation D = Do + X A for the quartiles in the LCRS, ACDM2 halos, ACDMl mock, and OCDM 
mock samples, and Poisson samples. 



samples lower quartile median upper quartile 



LCRS data samplings 


6.7± 1.1 ^ 


- (1.5 ± 0.1) 


X 


A 


14.9 ±2.0- 


f (2.1 ±0.2) X A 


26.0 ± 


3.0- 


ir (2.4 + 0.3) X 


A 


LCRS data for rotated grid 


8.5 ± 1.0 ^ 


- (1.4 ±0.1) 


X 


A 


18.1 ± 1.7- 


f (1.8 ±0.1) X A 


25.5 ± 


2.5 - 


^ (2.4 ±0.2) X 


A 


LCRS data for diluted samples 


7.4 ± 0.7 ^ 


- (1.5 ± 0.1) 


X 


A 


16.5 ±1.3 - 


f (2.0 ±0.1) X A 


28.4 ± 


1.9- 


^ (2.2 ±0.1) X 


A 


Base void sizes (Miiller et al.) 


5.7 ± 1.6 ^ 


- (0.9 ±0.1) 


X 


A 


11.8 ±2.9- 


f (1.1 ±0.2) X A 


16.8 ± 


2.9+ (1.5 ±0.2) X 


A 


ACDM2 FOF-halos 


7.5 ± 1.4 ^ 


- (1.5 ± 0.1) 


X 


A 


9.6 ± 1.2 4 


- (2.0 ±0.1) X A 


20.4 ± 


2.0 + (3.0 ±0.1) X 


A 


ACDMl mocks 


2.0 ± 0.9 ^ 


- (1.8 ± 0.1) 


X 


A 


6.8 ±2.2 4 


- (2.9 ±0.2) X A 


19.7 ± 


1.9 - 


^ (2.9 ±0.2) X 


A 


OCDM mocks 


4.3 ± 0.6 ^ 


- (1.7 ±0.1) 


X 


A 


11.2 ±0.9 - 


f (2.3 ±0.1) X A 


16.3 ± 


1.5 - 


^ (3.2 ±0.1) X 


A 


Poisson 


1.0 ± 0.6 ^ 


- (1.9 ±0.1) 


X 


A 


0.0 ±0.9 4 


- (3.0 ±0.1) X A 


1.1 ± 


1.0 4 


(3.8 ±0.1) X 


A 



we could not answer the question for the relation of the void 
statistics in the narrow LCRS slices to the void statistics in 
full 3D galaxy samples. Both questions are addressed in the 
present paper using a set of realistic simulations with dif- 
ferent procedures for galaxy identification in high resolution 
dark-matter simulations with different resolutions and box 
sizes. 

Here we are not basically interested in discriminat- 
ing between different cosmological models, therefore we re- 
stricted us mainly to the standard ACDM model with a den- 
sity parameter Qm ~ 0.3 and a Hubble constant h = 0.7, and 
we took as a comparison model only an open model OCDM 
with Qm = 0.5 and h = 0.6. These parameters correspond 
to the models discussed in our earlier paper, and therefore, 
they are well suited to answer the mentioned questions. In 
Table ^ we give the simulation details, in particular the nor- 
malization by the linear mass variance on a 8 /i~^Mpc scale, 
(jg, the box size, and the mass and spatial resolution. 

We perform two low resolution particle-mesh (PM) sim- 
ulations with 256^ particles in large boxes of (500/i"^Mpc)^ 
size using 512^ grid cells, they are denoted as ACDMl and 
OCDM. These simulations are siin ilar to them used previ- 
ously by ( [Doroshkevich et al. 1999| ) and Miiller et al. (2000). 
With the large volume, they are well suited for discussing 
the relation of voids in the effective 2D slices of the LCRS 
to the true 3D void statistics. In these models we identified 
galaxies with single particles employing a nonlinear bias pre- 
scription detailed in Miiller et al. (2000). In short, it employs 
the local density constructed from the neighboring particles 
and uses a threshold bias prescription for suppressing galaxy 
formation in low density regions, and a non- linear bias for 
particles in high density regions to model the merging and 
higher dark matter to galaxy ratio in high density regions as 
caused by the continuous merging in the hierarchical clus- 
tering process, cp. Cole et al. (1998). 

Further we employ medium resolution particle-particle 
particle-mesh (P3M) simulations with the adaptive mesh re- 
finement code of Couchmann (1991) using also 256^ parti- 
cles. These simulations, denoted ACDMl and ACDM2, use 
boxes of (280/i"^Mpc)^ and (200/i^^Mpc)^, respectively, to 
test for the mass resolution dependence and the influence 
of the upper box size on the void statistics. The mass reso- 
lutions are 10^^ and 4 x 10^° h~^MQ, respectively, and we 
employ a spatial resolution of 50/i~^kpc and 40/i~^kpc in co- 
moving coordinates, cp. Table This means we are able to 
resolve galaxy sized halos with 10-30 particles, but we also 



get larger group and cluster sized halos that enclose systems 
of galaxy halos being overmerged in the medium resolution 
simulations. For the question of voids, the distribution of 
galaxies within groups and clusters is of minor importance. 
Therefore, we can identify field galaxy halos with systems 
of about 10 particles and we study the void distributions 
within these systems. 

For the higher resolution simulation ACDM2 and 
ACDM3, we use a standard FOF halo finder, with a dimen- 
sionless linking length (in units of the mean particles separa- 
tion) fe = 0.12 (ACDM2) and fe = 0.17 (ACDM3). The rela- 
tively small linking lengths effectively collect galaxy like ha- 
los also in high density regions of the simulation as in groups 
and clusters which is essential to get a realistic indication of 
void boundaries in higher density regions of the cosmic web. 
We impose a mass boundary of 13 DM-particles for ACDM2 
and 8 dark matter particles for ACDM3, but for many anal- 
yses we further restrict the minimum of dark matter par- 
ticles per halo. After identifying the halos for ACDM2, we 
check the virial theorem of the FOF- groups, and we split 
off groups with high kinetic energy T as compared to the 
potential energy \W\, T > \W\, to separate interacting or 
merging groups. To this aim we employ a 6D-group finder 
with a linking measure Ar^/r^ + Av^ /a'^ where r and 
are the effective (half-mass) radius and the velocity disper- 
sion of the unbound groups. This concerns about 10% of the 
halos, with many halos having only a small number of DM 
particles. The splitting off of unbound halos adds about 2% 
to the equilibrium halos due to a fixed mass boundary for 
the final halos. Thereby we get about 53000 halos resam- 
pling possible galaxies in ACDM2 corresponding to a mean 
separation of A = 7.45/i~^Mpc. 

For ACDM3, we select halos at redshifts z = 0, 1, 2, 3 to 
check for the evolution of the void sizes. We take all halos 
selected with a linking length b — 0.17, finding about 160 
000 objects with a mass limit of 8 dark matter particles, 
i.e. with a mean interparticle separation of A = 3.7/i~^Mpc. 
At 2: = 3 this grows to A = 4.4/i~^Mpc due to the smaller 
number of halos at this redshift. We also employ a minimal 
circular velocity Vdrc > 200 km/s for identifying normal 
mass galaxy halos which leads to about 23000 objects with a 
mean separation of A = 7/i~^Mpc at z = 0. The mean sepa- 
ration between halos grows at higher redshift, corresponding 
to A = 7.8/i"^Mpc at z = 1, A = QAh'^Mpc at z = 2, and 
A — 13ft^^Mpc at 2 = 3 due to the mass increase in single 
halos by merger precesses and contirmous mass accretion. 
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Figure 3. Scaling relation for median and quartile values of the 
void sizes of mock galaxy samples of the OCDM simulation in 
a 500/i~^Mpc box (filled and open circles and dotted lines) and 
halos in the ACDM2 simulation using a 280/i~^Mpc box (stars, 
crosses, and dashed lines), both in a geometry similar to the 
volume-limited sub-samples of LCRS (solid lines with error bars) . 
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Figure 4. Comparision of the scaling relation in 2D and 3D for 
mock galaxy samples of the OCDM simulation in a 500/i-lMpc 
box. The 2D results are based on a geometry similar to LCRS 
volume-limited samples (small circles and dashed lines). The 3D 
analysis leads to smaller values of the medium and quartile void 
sizes (large squares and solid lines), but the slope of the scaling 
relation is similar. 



As shown in Muller et al. (2000) for the nonlinear bi- 
ased mock catalogues, the halos reproduce well the two-point 
correla tion function of gala xies in the LCRS in the redshift 
space (Tucker et al. 1997). Therefore they are reasonable 
test objects tor the void statistics in the underlying CDM 
models. 

The key questions of the simulations is whether they 
can reproduce the void size distribution, and in particular, 
the slope of the scaling relation for randomly diluted model 
galaxy samples. This characteristic measures the clustering 
properties in different parts of the universe observed by the 
LCRS. As it was shown in the previous paper, this slope is 
smaller for the LCRS samples than for the mock samples of 
N-body simulations, i.e. the median and upper quartile void 
sizes determined in the dense parts of the LCRS are larger. 
This was shown using mock samples having the same geom- 
etry as the LCRS samples. Especially in the well sampled 
parts the mock samples are not able to reproduce large voids. 
We complement this analysis by studying friend-of-friend 
dark matter halos of the ACDM2 simulation (cp. Table |^. 
In Fig. ^ we show that the scaling relation of voids in the 
OCDM mock samples (dotted lines) has almost the same 
steep slope as the scaling relation of halos in the ACDM2 
simulation (dashed line, cp. also Table |l]). In particular, it is 
significantly steeper than that from the LCRS (solid lines). 
The analysis of the halo samples was performed by con- 
structing mock samples in a slice geometry as in the LCRS 
data, and we did a 2D void search as described in the last 
section. To resume, we cannot reproduce the typical voids 
in well sampled parts of the LCRS as shown on the left of 
Fig-i 

The values found for the slopes of the scaling relations 
for the quartile and median values are given in Table 0. It 



shows that the OCDM and ACDMl mock samples and the 
ACDM2 FOF halos, that were constructed in different sim- 
ulations, have similar scaling relations, in particular similar 
slopes. In all samples, we notice a statistically significant 
discrepancy to the LCRS. 



5 VOIDS IN 3D AT REDSHIFT ZERO 

While in Section ^ mock and halo samples with a slice ge- 
ometry similar to the LCRS were analyzed, now we pro- 
ceed to 3-dimensional samples obtained from simulations. 
As a typical example results from the OCDM simulation in 
a 500/i^^Mpc box are shown in Fig. ^. In 10 different cuts 
through the box, diluted mock samples with LCRS geom- 
etry were constructed. Their scaling with the mean sepa- 
ration is indicated by the dashed lines. Additionally the 3- 
dimensional algorithm is applied to the complete box. The 
full sample was also diluted in this case to investigate the be- 
havior with lower sample density. Not only do the voids obey 
the scaling relation, but also the scatter is typically much 
lower. The scaling relations show the same slope, while there 
is an overall offset in the size. We regard this difference as a 
feature of the two dimensional slices. By observing only thin 
sheets in the underlying 3-dimensional galaxy sample, as it 
is the case with the LCRS slices, the effective size of voids 
increases due to missing galaxies in the bounded region of 
the slices. The increase in the absolute value of the scaling 
relation obviously depends on the mean depth of the slice 
in comparison to the mean void parameters, and we cannot 
connect it with any fundamental property of the galaxy dis- 
tribution. The observed void sizes become on average 20 - 
30% smaller in 3D, but the slope of the scaling relation of 
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Table 2. Cosmological A'-body simulations 



model 




h 


0-8 


box 


size 


mass resolution 


spatial resolution 


ACDMl 


0.3 


0.7 


0.87 


500/1- 


iMpc 


6 X 10" /i-IJVIq 


1.0 /i-^Mpc 


ACDM2 


0.3 


0.7 


0.91 


280/1" 


iMpc 


10" /i-iM0 


0.05 /t-^Mpc 


ACDM3 


0.3 


0.7 


0.87 


200/1- 


iMpc 


4 X lOio /i-iM0 


0.04 /i-^Mpc 


OCDM 


0.5 


0.6 


0.80 


500/1" 


iMpc 


10^2 /i-iM0 


1.0 /i-^Mpc 



void sizes with the mean galaxy separation remains constant 
and represents an important characteristics of the hierarchy 
of the galaxy distribution in the cosmic web. 

An important question to be addressed while study- 
ing void sizes is the dependence on the mass of the objects 
which represent the structure under analysis. From elemen- 
tary considerations typical voids expected within the distri- 
bution of dwarf galaxies as an example must be different in 
size from the voids present in the large scale structure built 
by galaxy clusters, as an extreme opposite example. This is 
firstly due to the different number densities of these objects, 
and secondly because smaller objects are distributed more 
smoothly on average than large objects. Both of these prop- 
erties affect the sizes of the voids, and they can be quantified 
by the scaling relation. 

In Miiller et al. (2000) we studied a possible dependence 
of median void sizes on the absolute magnitude. The magni- 
tude range of galaxies in the LCRS varied between -20 and 
-21. This is a small range and as such we could not find 
any significant change in the sizes of voids. Here we study 
the dependence of void sizes on the FOF halo mass. Fig. H 
shows that the cumulative void size distribution bends to- 
wards larger void sizes with growing mass of the halos. A 
ACDM simulation in a 280/i~^Mpc box with a mass resolu- 
tion of 1O"M0 was used. The lower cutoff in the mass range 
was set to 1.3 x 1O'^^M0 for the most complete sample and 
to 1.5 X 10^^ M© for the one with the lowest space density of 
halos. In the first case the sample contains all objects with 
a mass comparable to a milky way galaxy and higher. In the 
latter case the objects that remain in the sample are com- 
parable to small galaxy groups. Because the number den- 
sity of the identified halos decreases with growing mass, the 
structure of the samples is dominated by halos with masses 
near the lower end. Within the mass range of halos under 
investigation the median value of the void sizes increases 
from 25/i~^Mpc to twice this value. Obviously the void size 
distribution looks remarkably similar, cp. also Miiller et al. 
(2000). 

As a comparison of the scaling relation from sparse sam- 
pling and increasing mass, we construct randomly diluted 
samples with average mean separations A similar to the halo 
samples with different masses. The diluted sets contain a 
mixture of high and low mass halos. Because of the higher 
abundance of the low mass halos, they dominate the diluted 
samples. Althouth dilution increases the void sizes, the voids 
in the diluted mixed samples are significantly smaller. As 
it is shown in Fig. ^ the difference grows to 20% for the 
sample with the largest halo masses. This effect is expected 
and easily explained by the picture of gravitational insta- 
bility. Larger mass halos are more clustered according to 
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Figure 5. The cumulative void size distributions in the ACDM2 
simulation in a 280/i~^Mpc box. We show the size distri- 
butions for dark matter halos with masses increasing from 
left to right, from 1.3 X lO^^ /i-^M© (solid line), 3.2 X 10^^ 
/i-iM0(dashed line), 6.0 X lO^^ /i-iM0(dashed dotted line), 
9.0 X 10^2 /i-lM0(dashed three-dotted line) to 1.5 X 10^^ 
/i~^M0 (dotted line). 



the higher amplitude of the 2-point correlation function of 
galaxy groups as compared to galaxies. Large and massive 
halos reside in the nodes of the large scale structure net- 
work, while the less massive ones are still streaming along 
pancakes and filaments towards these nodes. Their distribu- 
tion is more homogeneous and less clustered. Table ^ com- 
pares the parameters of the scaling relation for these two 
differently selected series of halo samples. A similar change 
of the void size was found in Little & Weinberg (1994), but 
they did not derive the scaling relation of the void sizes un- 
der random dilution studied here. There were found larger 
voids if a biasing prescription for galaxy formation was ap- 
plied. With the present approach we quantify this relation 
of object masses and void sizes. 

An additional study illustrates the relation of the void 
sizes to their environmental and central density. For this 
purpose a smoothed density field was constructed in the 
dark matter distribution of the ACDM3 simulation in a 
200/i~^Mpc box. In the first step the dark matter parti- 
cles were distributed onto a grid with 0.5 /i~^Mpc per cell, 
using the cloud-in-cell method. After smoothing the density 
field with a Gaussian function [a — \ /i~^Mpc) a threshold 
overdensity of &th ~ 5.0 was applied. Every grid cell with 
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Table 3. The scaling relation in 3D for the lower quartiles, medians and upper quartiles of void size distributions in randomly 
diluted OCDM and FOF-halo samples and in halo samples with increasing mass. 



samples lower quartile median upper quartile 

OCDM mocks, randomly diluted 3.3 ± 0.2 + (1.5 ± 0.02) X A 5.2 ± 0.1 + (2.5 ± 0.01) X A 9.8 ± 0.5 + (2.9 ± 0.06) X A 

ACDM halos, randomly diluted 3.7 ± 0.2 + (1.5 ± 0.02) X A 7.7 ± 0.4 + (2.4 ± 0.04) X A 13.5 ± 0.8 + (3.0 ± 0.07) X A 

ACDM halos, increasing mass 1.5 ± 0.3 + (1.8 ± 0.02) x A 4.1 ± 1.5 + (2.9 ± 0.13) x A 6.2 ± 0.8 + (4.1 ± 0.07) x A 




Figure 6. The median, upper and lower quartile values of the 
void size distributions of DM halos with growing mass as de- 
scribed in Fig. 1^ (squares and dotted lines), compared to diluted 
samples using the lowest mass threshold (diamonds and solid 
lines). 



Figure 7. Averaged density profiles of under dense regions 
in a ACDM2 simulation with 200/i^^Mpc box size. The aver- 
age is taken in size ranges 36/i~^Mpc - 50/i~^Mpc (solid line), 
24/i-iMpc - 36/i-lMpc (dashed line), im'^Mpc - 2ih-^Mpc 
(dashed-dotted line) and 10/i~^Mpc- 16/i~^Mpc (dashed-3dotted 
line). 



a density below this value was regarded as containing no 
significant structure. After finding the voids their density 
profiles were computed by averaging in concentric spheres 
around the center of the base voids. The density profiles as 
shown in Fig. R were obtained by averaging in size ranges be- 
tween (36-50)/i"^Mpc, (24-36)/i~^Mpc, (16-24)/i~^Mpc 
and (10 — 16)/i~^Mpc. To obtain the average between voids 
of different sizes we use a relative radial coordinate r/R, 
where r is the distance from the void center and R is the 
corresponding effective void radius defined above. For the 
study of the environment of the voids we extend the aver- 
aging to relative radial coordinates larger than unity. They 
are calculated up to a sphere with r = 2R, the last sphere 
contains a volume eight times larger than the void. 

Two effects are obvious in the profiles shown in Fig. ^ 
The average underdensity at the center of voids depends on 
the void size. On average the central density of small voids 
is higher than that of large voids. This correlation refiects 
the inversion of gravitational clustering. While high peaks 
in the density fluctuations evolve to the knots in the large 
scale structure, large voids are formed in regions with a pri- 
mordial large scale underdensity. Also the profiles show the 
tendency that the environment of large voids is underdense, 
while small voids will form both in underdense and over- 



dense regions. This is typical for the hierarchical clustering 
scenario. 



6 EVOLUTION OF VOIDS 

A study of the void evolution using an adhesion model of 
structure formation was carried by Sahni et al. (1994). There 
a correlation of the void size with the value of the primordial 
gravitational potential was suggested, and it was derived a 
redshift dependence of void sizes, D{z) = Do /{I + z)^^^. 
More recently, Friedmann and Piran (2000) proposed a 
model that combines the growth rate of negative density 
perturbations with a lower probability of galaxy formation 
in these underdense regions. According to this model a typ- 
ical void with a size of 20/i~^Mpc corresponds to a 3a per- 
turbation so that they occur very rarely. This simple model 
can explain relatively large voids in the case of a ACDM 
universe, but it still fails to reproduc e the occurrence o f as 
large voids as observed in the LCRS ( Miiller et al. 200C ), or 
in th e ORS and IRAS Galaxy catalogues, ( |E1-Ad fc Pirim 
200^). 



Here we make use of the ACDM3 simulation in a 
200/i~^Mpc box. In addition to the resulting dark matter 
distribution at present epoch, we stored the snapshots of 
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Figure 8. Cumulative size distribution of voids among darlc mat- 
ter halos in a ACDM simulation with 200/i~^Mpc size. The solid 
lines refer to z=0, dashed, dotted and dashed-dotted lines refer to 
z=l, z=2, and z=3, respectively. The distributions to the left are 
obtained in the case of complete halo distributions, and the oth- 
ers belong to the case with a lower cutoff in the circular velocity 
i^circ = 200 km/s. 



the evolution aX z — 1, z — 2 and 2 = 3. With the FOF- 
algorithm described above, gravitationally bound dark mat- 
ter halos of were identified at each time step to study the ef- 
fects of void size evolution. A minimum of 8 dark matter par- 
ticles was chosen as a requirement for the halo identification, 
but we keep in mind that these low mass halos are not well 
resolved and quite unstable to disruption in the tidal fields 
of neighboring halos. These halos have a minimum mass of 
3 X 10^^ Mq. Furthermore we construc ted halo s amples with 
a minimum circular velocity Vcirc ~ \/ GM/rvir = 200 km/s 
within the virial radius r^ir- 

As Fig. ^ shows there is little sign of evolution up to a 
redshift of z = 2 if we consider the complete samples. This 
means that void sizes grow in the same rate as the Ifubble 
expansion. The network of high and low density regions in 
the universe has been formed in an early epoch of the cos- 
mic evolution. But on the other hand a significant evolution 
of the void distribution is observed when we consider halos 
larger than a certain mass or circular velocity. This implies 
that after the early setup of the large scale network the mat- 
ter flows along the walls and filaments, continues to merge 
and forms larger objects with time. At the redshift 2 = 1, 
the evolution of the number density is already visible in the 
void distribution. At redshifts higher than unity, large ha- 
los are still in the stage of forming. Their number density is 
significantly lower, and the voids in the spatial halo distribu- 
tion are increasing in comoving coordinates. The formation 
of more massive objects leads to a subdivision of voids as 
we approach redshift z = 0. 

The evolution of the underdense regions is accompanied 
by a flow of dark matter from low dense to high dense re- 
gions. Here we study the evolution of the density proflles as 
described above. The resulting averaged density proflles are 
shown in Fig. ^at different redshifts. The density proflles de- 
velop distinctly for different sizes. The value of the average 



density in the central parts decreases linearly in all panels 
with 5 = 0.18 per redshift unit. The upper left panel of Fig. 
^ corresponds to the largest underdense regions. Different 
from the other size ranges, the average density in the envi- 
ronment of such voids decreases during the evolution from 
z = 3 to 2 = 0. 

In the upper right and lower left panel of Fig. ^, which 
show the proflles of median size underdense regions, the av- 
erage density of the void environment remains nearly con- 
stant and roughly at mean cosmic density. In the lower right 
panel, which corresponds to smaller voids, the average den- 
sity in the environment grows with redshift. Fig. ^ illustrates 
the flow of dark matter from low into high density regions 
mainly surrounding the small voids. A correlation between 
regions with positive initial gravitat ional potential and un- 
derdense regions was also found in ( Madsen et al. 1998 ), in 
good agreement with our interpretation. There it was also 
found a decreasing average density in the underdense re- 
gions. 



T DISCUSSION 

The present paper includes an extension of our previous 
analysis of the void size distribution in the 2D samples of the 
LCRS. By comparison of large 3D boxes with mock 2D slices 
in similar geometry as the LCRS we can extend our previ- 
ous results. Although there is more scatter and uncertainty 
in the statistical measures of the void sizes in 2D, there is 
evidence that these results are physically relevant for the 
complete spatial distribution of galaxies. From the simula- 
tions we flnd a strong hint that the observed scaling relation 
of effective void sizes with the mean galaxy separation in the 
LCRS samples has a counterpart in the 3D structure. The 
effective size was deflned both in 2D and 3D as the diameter 
of the equal sized circle or sphere as the voids. The two differ 
only by a nearly constant offset, i.e. the 3D void sizes that 
occupy the same volume fraction (or area for 2D) are about 
20% smaller. We test the robustness of the scaling relations 
of the LCRS by evaluating the variances that the differently 
selected homogeneous samples, the differently oriented grids, 
and the dilutions introduce into the results. Other possible 
sources of systematic errors, such as the influence of the se- 
lection function and geometry and boundary effects, were 
investigated by our previous paper. The uncertainty due to 
these effects was found to be less than 10%. 

The next part of our analysis is devoted to the ques- 
tion of how the void sizes in the spatial distribution of some 
objects depend on the type of objects. Here we study the 
dependence for halos of milky way mass and above. Below 
this mass range there might be some inconsistency of the 
galaxy population in voids betwe en the observa tional data 
and models of galaxy formation ( Peebles 2001 ). As we re- 
strict our samples to higher mass halos, it is naturally to 
expect that the identifled voids are larger, just because the 
number density of halos decreases and the mean halo sep- 
aration becomes larger. An important result obtained with 
our analysis is that larger voids sizes are found for halos with 
larger mass. This effect is shown by comparing the scaling 
relation for samples with increasing halo masses and diluted 
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Figure 9. The evolution of average density profiles for voids in the same ranges as in Figure |?|. Each panel shows the evolution in one 
size range, from top left the largest voids to lower right the smallest ones. The profile curves in each panel correspond to different epochs, 
z=0 (solid lines), z=l (dashed lines), z=2 (dotted lines) and z=3 (dashed-Sdotted lines) and illustrate the evolution in each size range. 



subsamples of the original sample with equal mean density. 
In the latter case the voids are smaller and the slope of 
the scaling relation is shallower because a mixture of halos 
with different masses is less clustered than the comparable 
sample of only high mass halos. Also in this respect the 
scaling relation of void sizes proves to be a good measure 
for the fundamental clustering properties of a sample. Our 
new statistics extends the void probability function that was 
shown to provide a generating function for the set of higher 
order correlation functions (White 1979, Sheth 1996). It is 
more sensitive to the large scales in the cosmic network of 
galaxy structures. 

Furthermore we discussed the void density profiles. 
These were estimated by averaging in spherical shells around 
the void centers. They show characteristics which are typ- 
ically difi'erent for large and small voids. The central and 
environmental density of large voids are significantly lower 
than small voids. The underdensity around a large void en- 
compasses a larger volume than the void itself, while the 
environment of small voids represents in the average an over- 



dense ridge. We also examined the evolution of void profiles 
up to a redshift z = 3. The central underdensity decreases 
linearly with S = 0.18 per redshift unit. Also the environ- 
mental underdensity of large voids becomes deeper during 
the evolution. For average sized voids the environment den- 
sity remains constant near the mean matter density, while 
the overdensity ridge of small voids grows. Regarding the 
size distribution and volume fraction of voids we find little 
sign of evolution in catalogues containing small dark matter 
halos. However if we restrict the samples to larger mass halos 
comparable to the milky way, there is a significant evolution 
in their number density and consequently in the void dis- 
tribution. As they continue forming they subdivide existing 
voids and shift the void size distribution towards lower sizes 
in comoving coordinates. 

The exploration of the void size evolution and of the 
void scaling relation in new redshift catalogues and in more 
realistic model galaxy distributions than considered here 
represents a challenging prospect for the future. It should 
allow a quantitative understanding of the large scale envi- 
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ronmental dependence of galaxy formation and its bias that 
should be basic for resolving the discrepancy of the scaling 
relation of the void size distribution function in the LCRS 
and in CDM models. In this way we can find an answer 
to the question, whether this discrepancy is a fundamental 
problem for the standard ACDM model or whether it can 
be explained by a more realistic description of the galaxy 
formation than it was possible by a nonlinear bias and by 
our friend-of-friends halo finder. 
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